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The derivation of Marko and Siggia's interpolation formula 
for the force-extension relation of the Worm- Like Chain Model 
( C. Bustamante, J. F. Marko, E. D. Siggia, and S. Smith, Sci- 
ence 265, 1599 (1994); J. F. Marko and E. D. Siggia, Macro- 
molecules 28, 8759 (1995)) is retraced. Isotropy of space, re- 
sulting in rotational invariance of the free energy, is invoked 
together with analyticity. A new interpolation formula results 
for the force-extension relationship. It is as simple as the old 
one, but twice as close to the exact force-extension relation- 
ship. Application of the same reasoning to the second-order 
perturbative result obtained at strong force (J. D. Moroz and 
P. Nelson, Proc. Natl. Acad. Sci. USA 94, 14418 (1997)) re- 
sults in yet a new interpolation formula, good to 1% at all 
forces. 



Introduction. The worm-like chain (WLC) model 
is the quintessential model of entropic elasticity 
from a flexible, but unstretchable fiber, string, or thread. 
The model is conceptually simple, mathematically mini- 
malist, and widely used to interpret experiments that in- 
volve pulling at strings at finite temperature. Thus recent 
single-molecule experiments in biological physics cause 
new interest in this old model from polymer physics. 
It was successfully employed to model the experimen- 
tal force-extension relationship of double-stranded DNA 
some years ago. Also, the force-extension relations 
for the giant muscle protein titin ||,[l0]], the polysaccha- 
ride dextran , and single molecules of xanthan Jl^] 
were explained with the WLC model. A refined analy- 
sis of DNA's force-extension relationship in terms of the 
WLC model was recently presented in The relax- 
ation dynamics of extended DNA molecules, measured 
with millisecond resolution and femtonewton force spec- 
troscopy, was interpreted using the WLC model in p4[ . 
In the WLC model was used to interpret stretch- 
ing of single collapsed DNA molecules. Furthermore, the 
WLC model was extended with stretch |8|,[l3| and twist 
1010 — 01 ancl references therein] to model also these 
property of double stranded DNA. 

The WLC model. The WLC- models describes a 
string of vanishing cross section, unstretchable, but flex- 
ible. As it cannot stretch, the string can only bend, and 
it resists even that with a force (per unit length) propor- 
tional to its curvature. The constant of proportionality, 



A, is called the bending modules, and has dimension en- 
ergy per unit length of string, per unit curvature squared. 
Thus the bending energy of the string is 



r ^ A f L ° , f dt 
Ei] = - ds[ — 



(s) 



(1) 



where Lq is the length of the string; s parametrizes points 
on the string by their distance from one end, as measured 
along the string, s £ [0, L }; t(s) is the tangent vector to 
the string at its point at s. This energy is evidently 
minimal for a straight string, since a straight string has 
constant tangent vector. 

When the string is submerged in a heat bath, a bend- 
ing energy of order /cbT is available to each of its de- 
grees of freedom. The string consequently bends in a 
random manner obeying Boltzmann statistics, and an 
attempt to pull apart its ends is resisted with a force F 
which depends on the string's end-to-end separation R, 
F = F(R). Because of the isotropy of space and the 
rotational invariance of the bending energy in Eq. (Q), 
F(R) is anti-parallel to R, and |F| depends only on 
R = \R\. Because the string is unstretchable, the largest 
possible end-to-end separation is R = Lq. So a natu- 
ral dimensionless measure of the end-to-end separation 
is r = R/Lq £ [0,1]. 

Despite its simplicity, the WLC model cannot be solved 
analytically in general. When its ends are left free, a 
small calculation based on the Boltzmann weight fac- 
tor exp(— E\t\/k-QT), results in the correlation function 
(t(si) ■ f(s2)) = exp(— \s\ — S2\/L p ), where the persis- 
tence length L p is inversely proportional to the temper- 
ature: L p = A/k^T J5j. In the limit L p /Lq — > oo where 
the string is much longer than its persistence length, 
L v F{r)/k^T is a dimensionless function of the dimen- 
sionless variable r, only. But not even in this convenient 
limit is an exact analytical solution possible. A numer- 
ical solution for the force-extension relationship is not 
difficult to obtain [[7p^] , and is given in a useful form in 
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The Marko-Siggia interpolation formula. It is some- 
times convenient, however, to have a simple analytical 
expression for the force-extension relationship, even if 
only an approximate one. Marko and Siggia presented 
such a relationship in 
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This formula was derived by calculating the force- 
extension relationship analytically to leading order in the 
limit of very large force, where the string is nearly fully 
stretched, r sa 1 @. That yielded 



kn~f 



1 



4(1 -r) 2 



unknown subdominant terms. 
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An attempt to use the known part of this result for all 
forces/all r 6 [0, 1] fails at small forces/small r, where a 
calculation Ml shows that 
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The explicit term in Eq. (J3J) does not satisfy Eq. (Q). 
The last two terms on the right-hand-side of Eq. (ph were 
added to ensure Eq. (||) is satisfied. We note that this 
procedure does not compromise the validity of the result, 
Eq. at r — > 1, because the two terms added remain 
finite in that limit. The result is accurate to 17% when 
at its worst — see Fig. 1 — and is by construction asymp- 
totically correct for r — > and for r — ► 1. 

Another interpolation formula. Now consider the ex- 
act force-extension relationship. Suppose we could calcu- 
late the string's free energy, T, analytically at given tem- 
perature and end-to-end separation R. Then we could 
calculate its force-extension relationship as 
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We cannot do this, but we know that the free energy is 
independent of the direction of R. It depends only on 
R = \R\. Furthermore, we expect the free energy to be 
an analytical function of R for \R\ < Lo, hence analytical 
in R = 0. Arguments for analyticity may be given |22| , or 
one may regard analyticity as a conjecture or postulate. 
Or one may simply disregard the issue; the interpolation 
formulas given below have the precision demonstrated in 
Fig. 1 no matter how we arrive at the formulas. Faulty 
logic works fine here. 

Analyticity of F(R) in R = implies that F(R) has 
a Taylor series expansion in powers of R's components. 
Rotational invariance consequently implies that R's com- 
ponents only occur in the combination R 2 in this Taylor 
series. Consequently, T is an analytical function of R 2 , 
T = F(R 2 ). Thus T is an even analytical function of 
R. From Eq. (||) then follows that F(R) must be an odd 
analytical function of R. 

Marko and Siggia's interpolation formula, Eq. (j^), is 
not an odd analytical function of r. But if we retrace 
its derivation from Eq. (Bl), we note that this expression 



already contains what it takes to mend it: we extrapolate 
Eq. (||) from the limit r — > 1 to lower values of r in an odd 
manner by realizing that the factor 1 in the numerator 
really is r, while the factor 4 in the denominator is (1 + 
r) 2 . Thus we arrive at a new analytical interpolation 
formula for the force-extension relationship of the WLC 
model: 
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The last term on the right-hand-side has been added to 
ensure that Eq. (||) is satisfied, entirely in the spirit of 
Marko and Siggia's derivation of their formula. 

i*8%(r) should be a better approximation than -FMs( r )> 
because it respects rotational symmetry and correctly 
captures all odd terms in the Taylor series for the ex- 
act result. They vanish in the exact result, and they 
vanish in F s %(r). Figure 1 illustrates this improvement: 
F s %(r) is much closer to the exact result than -Fms(?") 
is. Its differs less than 8%, at most, while -FmsM differs 
up to 17%. At low values of r, F 8 %(r) represents an im- 
provement by much more than a factor two. For larger 
values of r, the improvement is approximately a factor 
two. Note that F 8 % achieves this doubled precision with 
half as much "patching" at r ~ 0: only one correcting 
term was added. 

Yet another interpolation formula. Moroz and Nel- 
son have calculated the first correction term to Eq. (|^) 
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This result also is not an odd function of r. But again this 
is easily remedied. The first term on the right-hand-side 
was treated above, and the next two terms are treated 
similarly, yielding 
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This last expression is an odd analytical function of r, 
but it does not satisfy Eq. ([!]) , and we cannot mend that 
simply by adding terms which remain finite for r = 1, 
as Marko and Siggia did. But we can proceed entirely in 
their spirit, and add terms of same order as the neglected 
terms, in the present case O (r(l — r 2 )). Doing that, we 
arrive at 
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Figure 1 shows that F^ir) reproduces F exact to within 
1.5%. 

Discussion. It is clear from the procedure we have 
used that one may continue it systematically by calcu- 
lating more terms in the two series for -FexactWs asymp- 
totic behavior at r = and r = 1, respectively. This ex- 
act asymptotic information can then be incorporated in 
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FIG. 1. F M s/foxact (upper 
curve), and F^/Fe^t (bottom 
maximum of JMs/^exact is 1.17, 
maximum of i ? g %/J ? oxac t is 1.08, 
maximum of Fi%/F exac t is 1.015 
ratios plotted here are identical 
persistence length suggested 
comparisons. 
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curve), F 8% /Fex a ct (middle 
curve) plotted against r. The 
and occurs at r = 0.56. The 
and occurs at r = 0.64. The 
and occurs at r = 0.66. The 
to the inverse of the effective 
pi Eq. (14)] for this kind of 



an increasingly complex result, by including an increas- 
ing number of terms of the general form r 2n+1 (l — r 2 ) rn , 
with suitable coefficients and exponents n and m. In 
view of the accuracy already achieved with Eq. (^) , this 
is hardly worthwhile for most purposes. 

For one purpose, however, it looks promising: as a 
way to present a high-precision analytical interpolation 
formula meant for numerical evaluation ]25| ]. At weak 
and strong force (small and large end-to-end separations) 
where numerical methods typically fail unless special care 
is taken, such a result is exact to a chosen order in pertur- 
bation theory. And chosen properly, this order renders 
the result uniformly good to a desired precision, for all 
forces/end-to-end separations. Bouchiet et al.'s numeri- 
cal interpolation formula (l]]2(J] is of this nature, being 
exact to leading order at small and large force, and uni- 
formly good to 1% p(| . 

The improved force-extension formulas presented here 
remain valid when the WLC model is extended to de- 
scribe a somewhat stretchable string as done in Q . This 
because the extension consists in allowing Lq to stretch in 
a simple manner depending on R, while F(r) = F(R/Lq) 
is left unchanged. 

It may be of interest to apply the approach used here to 
the WLC model extended with twist. This model is ana- 



lyzed with strong- force perturbation theory in ]1J,|18|,|19| . 
So the results obtained there might have their range of va- 
lidity for a given precision extended down to lower force. 
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